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Abstract. The inverse Kasteleyn matrix of a bipartite graph holds much information 
about the perfect matchings of the system such as local statistics which can be used to 
compute local and global asymptotics. In this paper, we consider three different weight- 
ings of domino tilings of the Aztec diamond and show using recurrence relations, we can 
compute the inverse Kasteleyn matrix. These weights are the one-periodic weighting where 
the horizontal edges have one weight and the vertical edges have another weight, the q vo1 
weighting which corresponds to multiplying the product of tile weights by q if we add a 
'box' to the height function and the two-periodic weighting which exhibits a flat region with 
defects in the center. 



1. Introduction 

1.1. Terminology. Domino tilings of bounded lattice regions have been extensively re- 
searched during the last twenty years. These tilings are the same as perfect matchings of a 
bounded portion G of the dual square lattice, in the following way: a matched edge corre- 
sponds to a domino; the fact that the dominos do not overlap means that no two matched 
edges share a vertex, and the fact that the dominos cover the entire region means that each 
vertex in the region is covered by a matched edge. In the statistical mechanics literature, 
one speaks of dimer covers rather than perfect matchings, and dimers rather than matched 
edges. 

The most well-studied example of such a model is domino tilings of the Aztec diamond 
which was introduced in [11] . Here, one tiles the region {(x,y) : \x\ + \y\ < n + 1} with 
2 by 1 rectangles where n is the size of the Aztec diamond. There are other examples of 
the theory, but they involve replacing the graph G with a different one, such as the regular 
square-octagon lattice (giving the so-called diabolo tilings) or the hexagonal mesh (giving 
lozenge tilings). 

By giving each edge a multiplicative weight, we can consider random dimer coverings: the 
probability of each covering is proportional to the product of the edge weights of the dimer 
covering. The corresponding discrete probability space is called the dimer model. If the 
graph G is bipartite (as it shall be for the rest of this paper) then each dimer covering can 
be encoded by a three dimensional discrete surface, where the third co-ordinate is derived 
from the specific dimer covering and is called the height function [30J . For random tilings, 
[HI E] showed that with probability tending to one, the height function of a randomly tiled 
large bounded region tends to a deterministic limit shape. This shape is not smooth over 
the entire region: typically, there are macroscopic regions wherein the tiling is "frozen" (i.e. 
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exhibits deterministic correlations), called facets; as such the measure is often said to be in a 
solid state here [22]. Outside of the facets, the correlations between pairs of dimers are only 
mesoscopic, tending to zero as the dimers move farther apart. If the decay is polynomial, 
the measure is said to be liquid; if it is exponential it is said to be gaseous [22J. Not all 
tilings possess a gaseous region; however, all but the most degenerate have liquid regions. 
The limiting height function is smooth in these regions. Figure [T] shows two random tilings 
of relatively large Aztec diamonds. 

1.2. Local asymptotics for nice regions. For particular bounded regions, one approach 
to study these models uses an interlaced particle system which can be derived from the 
underlying tiling |15|,l3]. Using the Lindstrom-Gessel-Viennot theorem [161 129] combined with 
the Eynard-Mehta theorem (e.g. see [1]), it is often possible to find the correlation kernel for a 
determinantal process and compute finer statistics for the model [T5| 3J. Using these statistics, 
one can study the fluctuations between the interface of the solid and liquid regions when 
the system size gets large. Amazingly, these fluctuations have the same distributions arising 
from the study of eigenvalues of random matrices (see for example (15 |, I1T| . fTH] . Furthermore, 
one can even change the boundary conditions of the underlying tiling problem to find more 
degenerate kernels which also appear in the random matrix literature, for example, see [T]. 

For bipartite graphs, the Kasteleyn matrix is a signed weighted adjacency matrix indexed 
by the white and black vertices of the graph [18] . The inverse of the Kasteleyn matrix, known 
as the inverse Kasteleyn matrix, for bipartite graphs provides much information about the 
model - by |19j the edges form a determinantal process with the kernel given by the inverse 
Kasteleyn matrix. Hence, by knowing the inverse Kasteleyn matrix for a bipartite graph 
one can compute all finite, local and global asymptotics of the edges in the dimer model. 
For lozenge tilings, the interlaced particle system kernel can be used to compute the inverse 
Kasteleyn matrix as the particle system kernel and the inverse Kasteleyn matrix are in 
bijection [26J. However, for domino tilings on the Aztec diamond, the most natural kernel 
from the interlaced particle system contains different information to the inverse Kasteleyn 
matrix does. Although, one can derive the particle system correlation kernel from the inverse 
Kasteleyn matrix, the particle system correlation kernel gives a better description of the 
interface between the unfrozen and frozen regions, see [7]. By knowing both the interlaced 
particle system and the inverse Kasteleyn matrix, we believe that one understands the full 
asymptotic picture of the system. 

1.3. Purpose. The aim of this paper to highlight an elementary procedure which allows 
one to compute the correlation kernel of the determinantal process associated to the edges 
of tilings of Aztec diamonds. For pedagogical reasons, we do this first for the most well- 
understood instance of the dimer model (thereby recovering the work of |14j). followed by 
two, substantially more complicated new settings: the so-called g vo1 weighting (Section [4J 
and the two-periodic weighting [5] which includes as a special case the uniform measure on 
diabolo tilings on a fortress graph [27J. The two-periodic weighting has a 2 by 2 fundamental 
domain which is defined in [22] . 

For the g vo1 weighting, large random tilings of the Aztec diamond possess a limit shape 
exists when q — > 1~ as the system size tends to infinity which exists using the results from 
|21j . However, when q — )■ l - and the a — > 00, the results from |21j will no longer apply but 
simulations seem to suggest that there may be a limit shape and possibly interesting local 
and global asymptotic behavior. Figure ^1 shows relatively large tilings with q vo1 weighting. 
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Large random tilings of the two-periodic weighting of the Aztec diamond feature all three 
phases where the gaseous region described by an 'octic' curve. Explicit formulas for this curve 
are given derived in |27j and [21] using general machinery. Figure [3] shows a relatively large 
random tiling of a two-periodic weighting of the Aztec diamond. By having an expression 
for the inverse Kasteleyn matrix for the two-periodic weighted Aztec diamond, it may be 
possible to study this model on all the phase interfaces. Computing the interlaced particle 
system kernel for the two-periodic weighting of the Aztec diamond using the Lindstrom- 
Gessel-Viennot theorem [16} 129] combined with the Eynard-Mehta theorem (e.g. see [I]) 
seems somewhat complicated - one can either proceed by inverting either a block LGV 
matrix or block Toeplitz matrix. 

1.4. Explicit inversion of Kasteleyn matrices. We turn now to a discussion of our 
methods. It is possible to derive the inverse Kasteleyn matrix for domino tilings of the 
Aztec diamond with weight 1 for horizontal tiles and weight a for vertical tiles (one-periodic 
weighting). We also show that it is possible to extend the method to the (/-analog of domino 
tilings of the Aztec diamond (q vo1 weighting) . The height function of a domino tiling can be 
viewed as the surface of a pile of Levitov blocks [231 124] : one can modify the edge weights so 
that the weight of each tiling is proportional to (^{Levitov blocks}. Finally, we show that 
it is possible to compute the inverse Kasteleyn matrix with a two-periodic weighting of the 
Aztec diamond and do so for an Aztec diamond of size Am for m £ N. 

It is often quite difficult to invert a matrix whose entries have parameters; indeed, the 
typical methods in the literature involve first orthoganalizing, so that the matrix to be 
inverted is diagonal [5]. However, in the particular case of the Aztec diamond, one can make 
some progress by computing a multivariate generating function for the entries of K . This 
is possible essentially using graphical transformations similar to those used in |11] which is 
the same procedure used in the generalized domino shuffle [27] to compute the weights in 
the transitional steps of the shuffling algorithm. 

We use three relations among the entries of K^ 1 . The first two are the matrix identities 
K.K^ 1 = I, K^ 1 .K = I; the third is a recurrence relation in n, the order of the Aztec 
diamond, generated from the generalized domino shuffling algorithm which determines only 
those elements of K -1 which correspond to two dimers on the boundary of the region. We 
solve this recurrence by computing the ordinary power-series generating function for its 
coefficients, which we call the boundary generating function. 

Having done this, the entries of K" 1 are slightly overdetermined. As such, we treat the 
equations K.K^ 1 = I, K -1 .K = I as linear recurrence relations, for which the boundary 
generating function serves as a boundary condition. 

1.5. Details. The heuristic in the preceding sections drastically oversimplifies the intricacy 
of the combinatorics involved in these calculations, which consume the bulk of this paper. 
We also found that the computation of the boundary generating function of the two-periodic 
weighting was intractible without the use of computer algebra. Here is a description of the 
calculation in somewhat greater detail. The recurrences K.K -1 = K^ 1 .K = I allow us 
to write each entry of the inverse Kasteleyn matrix as a linear combination of entries of 
-K" _1 (w, b) where b and w are black and white vertices on the boundary. To find -fT -1 (w, b) 
where b and w are black and white vertices on the boundary, we first find Z(v,h)/Z where 
Z is the partition function (the sum of the weights of all tilings) and Z(v, b) is the partition 
function of the tiling with the vertices b and w removed. 
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As removing two vertices w and b from the boundary does not change the so-called Kaste- 
leyn orientation, we are able to recover K -1 (w, b) from Z(v,\))/Z. The boundary recurrence 
relation gives a relation for Z(w,b)/Z as a linear combination of Z(wi,b\)/Z for a smaller 
Aztec diamond where w, w\ and b, b\ are white and black vertices on the boundaries of 
their corresponding Aztec diamonds. To solve this recurrence, we use a generating func- 
tion approach and hence find -fT _1 (w, b) for b and w are black and white vertices on the 
boundary. However, for the q ml weighting, it is simpler to solve the boundary recurrence re- 
lation directly using double contour integral which we obtained from its boundary generating 
function. 

For the one-periodic and q weighting, we can use the three sets of recurrence relations 
mentioned above to find a formula for the inverse Kasteleyn matrix. For the one-periodic 
weighting, we write the formula as a generating function. For the g vo1 weighting, we obtain 
a double contour integral. Indeed, sending q — > 1, we find a double contour integral formula 
of the inverse Kasteleyn matrix in the one-periodic case, see [7] for a comparison. 

For the two-periodic weighting of the Aztec diamond, we must use a further two recurrence 
relations; K*K.K~ X = K*I and K~ l .K.K* = I.K* due to the complexity of the model; of 
course, these are completely equivalent and in principle carry no new information; however, 
they are simpler in a certain sense as they involve the discrete Laplacian. Additionally, the 
boundary recurrence relation for the two-periodic weighting of the Aztec diamond has order 4 
with the size of the Aztec diamond and is dependent on parity of the white and black vertices 
on the boundary. This leads to a matrix equation for the recurrence relations describing the 
boundary generating function. We only find the solution of the matrix recurrence relation 
when the size of the Aztec diamond is equal to 4m, though in principle the "only obstacle" 
to handling the other sizes of diamond is the difficulty of the computations. This leads a 
formula for the inverse Kasteleyn matrix when the size of the Aztec diamond is equal to 4m. 



1.6. Remarks about Asymptotics. As remarked above, this paper might provide a gate- 
way to computing fine asymptotics of the g vo1 and two-periodic weightings of the Aztec 
diamond. That is, one will hopefully be able to compute global correlations and local corre- 
lations at the phase boundaries using the formulas found in this paper. 

From simulations, the g™ 1 weighting seems to exhibit interesting limiting behavior when 
q — > 1~ and also when q — > 1~ and a — > oo as n — > oo simultaneously. We believe that a 
possible parametrization would be to set a = e C//2 and q = e~ c l n where n tends to infinity. 
Although it is possible to derive to the limit shape using [21] . we think our formula could 
be used to find the height fluctuations in the unfrozen region when c > 0. When c = 0, the 
model is equivalent to the one-periodic weighting of the Aztec diamond and so the height 
fluctuations are governed by the so-called Gaussian Free field (details of this process can 
be found in [28]) as shown in [7J. When c tends to infinity, we suspect that the unfrozen 
region is a flat square given by alternating (diagonal) columns of east and west dominos. 
That is, for a rescaled Aztec diamond with corners given by (0, 0), (0, 1), (1, 0) and (1, 1), the 
unfrozen region is given by (1/2, 0), (0, 1/2), (1/2, 1) and (1, 1/2). From initial computations, 
the asymptotic analysis to find these height fluctuations is encouraging and is current work 
in progress. 

It may also be possible to derive the q vo1 correlation (particle) kernel using established 
techniques (e.g. This works quite cleanly, since the process in question is Schur [2*5] , 

although, we have not tried this computation. 
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Figure 1. Random Simulations of the Aztec diamond of size 100 with one- 
periodic weight for a = 1 (left) and a = 1/2 (right) (see Section [3] for a 
description of edge weights) . 





For the two-periodic weighting, the process does not appear to be Schur. As mentioned 
above, under the right choice of parameters the model exhibits a third phase which has been 
named gaseous in which the height function correlations decay exponentially [22]. Other 
models that might possess similar phenomenon are the three periodic lozenge tiling in a 
large hexagon and the six vertex model with domain wall boundary conditions away from 
the so-called free fermion point |10j . Indeed, the six vertex model at its free fermion point 
with domain boundary conditions can be recovered from the one periodic Aztec diamond 

The main motivation behind this paper was to find the correlation kernel of the two- 
periodic weighting of the Aztec diamond, so that one can find the local correlations at the 
gaseous-liquid boundary. As mentioned above, this boundary represents the transition from 
the correlations of dominos having power law decay to the correlations of dominos having 
exponential decay. For tiling models, the solid-liquid boundary (with no cusps) has been 
well studied (see [26] for the most general results); the interlaced particle system associated 
with the tiling has fluctuations of size n 1 ' 3 and the distributions of particles are normally 
given by the so-called Airy process, a natural distribution originally formulated in Random 
matrix theory [2] . As far as we are aware, the gaseous-liquid boundary of a tiling model has 
not been previously studied in any such probabilistic model. 

The formula for the inverse Kasteleyn matrix for two-periodic weightings of the Aztec 
diamond, is given as a four variable generating function and it is not immediate how to 
analyze asymptotically. However, Kurt Johansson, using our formula, has been able to 
derive a double contour integral formula. From this double contour integral formula, it 
should be possible to use a saddle point analysis. Indeed, early computations using this 
approach show that the limiting octic curve can be recovered which agrees with limiting 
curve computed using the techniques from [21]. We believe that this approach will lead to 
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Figure 2. Random simulations of the Aztec diamond of size 200 with 
weighting (see Section [4]). The top picture has weights q = 0.99 and a 
The picture below has q = 0.98 and a = 10. 
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Figure 3. Random simulations of the Aztec diamond of size 200 with two- 
periodic weights (see Section [5] The first picture has a = 1/2 and 6=1 with 
8 colors. The second picture is the same tiling as the first but contains four 
colors to highlight the three phases for black and white printing. 
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Figure 4. The Kasteleyn co-ordinates for an Aztec diamond of order 3. 
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finding the gaseous-liquid fluctuations. This is current work in progress which will hopefully 
appear elsewhere [6]. 

1.7. Overview of the Paper. The paper is organized as follows: Section [2] we give some 
of the prerequisites and notation for understanding the proofs of the rest of the paper. 
In Section [3j we compute the generating function of the inverse Kasteleyn matrix for one- 
periodic Aztec diamonds which provides a blueprint for the computing. In Section |4j we give 
a contour integral formula for the inverse Kasteleyn matrix for q vo1 weighting of an Aztec 
diamond. Finally in Section [5] we derive the generating function for two-periodic weightings 
of the Aztec diamond for size 4m. 
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tion of computing K^ 1 for diabolo tilings which prompted this research and the very useful 
discussions which followed, Kurt Johansson for valuable discussions and improvements to 
this paper. We would also like to thank Alexei Borodin, Cedric Boutillier, Maurice Duits, 
Harald Helfgott, Richard Kenyon, Anthony Metcalfe and Andrei Okounkov for some useful 
insights. We would particularly like to thank MSRI Berkeley, where this work was initi- 
ated, and the Knut and Alice Wallenberg foundation grant KAW2010.0063, which partially 
supported the authors' work during its completion. 

2. Notation and Background Information 

Let W = {(xi, X2) : x\ mod 2 = 1, X2 mod 2 = 0, 1 < x\ < 2n — 1, < Xi < 2n} and 
let B = {(x\,X2) ■ x\ mod 2 = 0,X2 mod 2 = 1,0 < x\ < 2n, 1 < x<i < 2n — 1}. The 
set WU B denotes the vertex set of the dual graph of the Aztec diamond (rotated by —it/ 4 
and translated), where W denotes the set of white vertices and B denotes the set of black 
vertices. We call the above co-ordinate system (of the dual graph) of the Aztec Diamond, 
the Kasteleyn co-ordinates - see Figure [4] for an example. Let e\ = (1, 1) and e2 = (—1, 1). 

For a planar bipartite graph G = (V,E), the Kasteleyn orientation is a choice of edge 
weights so that the product of edge weights around each face is odd. In this paper, G 
is the graph formed from (the dual graph of) the Aztec diamond. We consider the edge 
weights given by positive numbers for edges parallel to e\ and positive numbers multiplied 
by i = for edges parallel to e2- As G is a bipartite graph, we consider the Kasteleyn 
matrix, whose rows are indexed by black vertices and columns indexed by white vertices, 
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with entries given by 

(2 1) K(b w) = { wt ^ iJ _1 for e = ( b > w ) and b - w = e i5 j e {1,2} 

^ ' ' ^ ' ' 1 otherwise 

where wt(e) is edge weight of edge e. As G is a bipartite finite graph, | det K\ is equal to 
the number of weighted dimer covers of G. This was first proved by Kasteleyn [18] but in a 
more general setting. The explicit formula for the inverse Kasteleyn matrix and correlation 
kernel of the dominos is given as follows: suppose that e\ = (pi, W\), . . . , e n = (b n , w n ) then 
the probability of seeing a perfect matching with the edges e\, . . . ,e n is given by [T!5] 

(2.2) det (KOh^K-^wM)^ . 

That is, the edges form a determinantal point process with the correlation kernel given 
by L(ei,ej) = K(bi,Wi)K~ 1 (uii,bj) [19]. In the case of the two-periodic weighting of the 
Aztec diamond, this point process is a block determinantal point process. With an explicit 



formula for the inverse Kasteleyn matrix, equation (2.2) means that we can compute any 
joint probability of a subset of edges of the matching which includes the probability of an 
edge. 

We now summarize the graph theory techniques used in this paper. 

Two dimer models are said to be gauge equivalent if their partition function differs by 
a constant multiple. The dimer model obtained from multiplying all the edge weights sur- 
rounding one specific vertex by the same constant is called a gauge transformation. Gauge 
transformed dimer models are gauge equivalent. 

Other than gauge transformations, we make use of three other graph transformations. All 
three of these alter the graph itself, but leave the partition function invariant up to a gauge 
transformation. 

(1) Suppose we have a large square with edge weights a, b, c and d where the labelling 
is clockwise, contained in some graph H. Suppose we deform this large square to a 
smaller square with edge weights A, B, C and D (same labelling convention as the 
large square) and also include an edge, with edge weight 1, between each vertex of 
the smaller square and its corresponding original vertex — see Figure [5] We call 
this new graph H*. We set A = cj (ac + bd), B = d/ (ac + bd), C = a/ (ac + bd) and 
D = bj (ac + bd) so that the local configurations and the weights of the matchings of 
H are preserved under the transformation to H*. For example, a dimer covering the 
edge with weight a and no dimer covering the edge with weight c in H corresponds 
to dimers covering the edge with weight C and the two diagonal edges incident to the 
edge with weight A in H* . This graphical transformation is called urban renewal [27] 
and we have that 

Zh = {ac + bd)Zn* 

where Zh and Zh* are the partition functions of H and H* respectively. 

(2) If a vertex, v, is incident to two edges, each having weight 1, we can contract the two 
incident edges and vertices of v, to v. The new edge set of v is the union of the edge 
set of the two contracted vertices omitting the two edges previously incident to v, 
all with the same edge weight before the contraction. This procedure is called edge 
contraction and has no effect on the partition function of the dimer covering. 

(3) If a vertex is incident to one edge, i.e. it is a pendant vertex, then this vertex, 
the incident edge and its incident vertex can be removed. This procedure is called 
removal of pendant edges. The partition function of the new graph formed by this 
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Figure 5. The Urban renewal step which maps the large square on the left 
the smaller square on the right and multiplies the partition function by ac+bd. 
The diagonal edges have weight 1 
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procedure is equal to the original partition function divided by the edge weight of 
the removed edge. 

The above three procedures can be used to compute the partition function [11], the edge 
probabilities [H [27] and the inverse Kasteleyn matrix for one-periodic, two-periodic and q vo1 
Aztec diamonds. For general weightings of the Aztec diamond and also other stepped square 
lattices, although it may be theoretically possible to find the inverse Kasteleyn matrix, we 
were not able to find suitable recurrence relations. 



In this section, we derive the inverse Kasteleyn matrix for uniform tilings as it provides 
a simplest example for our new method. The Inverse Kasteleyn matrix for the one-periodic 
weighting of the Aztec diamond was originally computed for uniform tilings in [T3]. In [7J, 
the formula was generalized to a biased one-periodic weighting using a one-to-two lifting 
from the interlaced particle system |15| . The new results here are: 

• the boundary generating function for the Aztec diamonds, defined below, and 

• the observation that the generating function for the entries of K -1 , for all Aztec 
diamonds together, is a rational function in six variables: two variables marking the 
row of K" 1 , two marking the column of K , one marking the size of the Aztec 
diamond, and the parameter a. 

For the setup, it is just as easy to consider the Kasteleyn matrix with horizontal tiles 
given weight 1, and vertical tiles given weight a (later, we will take a = 1 for simplicity). 
The Kasteleyn matrix for the one-periodic Aztec diamond is given by 



3. One periodic case 




1 if x — y = ±ei 
ai if x — y = ±e 2 
otherwise 



for x GB and y G W. 

Let f n {x) = (1 - x n )/{l - x) and let 





(wi,u>2) and b = (61,62) set 



(3.3) 



(3.2) 



F n °<°(w,b,a) = F n (wlb 2 2 ,a)w l b 2 , 
F°> 2 >,b,a) = F n (-l/ W l-bla)w\ n - l bl n b 2 i/a, 
F 2 "'°(w, b, a) = F n {-wl -1/61, ajw^b^i/a 
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and 

(3.4) F 2 "> 2 >,b,a) = F n {\lw\^lb\a)w^- X w^bf- X . 
Let 

(3.5) G n (w,b)= Yl ^K-\x,y) 

where w 1 = wf 1 w^ 2 and h y = b^ 1 6| 2 denote the generating function of the inverse Kasteleyn 
matrix of the one-periodic domino tilings of the Aztec diamond with Kasteleyn orientation 
given by multiplying all vertical edges (the vector e 2 is vertical) by i. 

Theorem 3.1. 



G„(w,b) 



w 1 w 2 b 2 w 2 f n+ ±(wlbj)f n (wjb 2 l ) 



( 3 - 6 ) + (1 + ±awi 



C(w 1 ,w 2 ) 

(1 + ia6 2 )F n °'° + b\ (b 2Wl f n {blwl) + (ia + b 2 2 )F^ 2n 
' C(«;i,«;2)C(6i,62) 



M 2, 2 blbl n+1 w lW rfn(blwj) + (1 + ia6 2 )F 2 "'° + b\{la + fe 2 )F 2 "' 2 " 

+ Ua + Wl)W2 C(w u w 2 )C( bl ,b 2 ) 

where C(ri, f2) = 1 + r 2 r| + ia(r 2 + r 2 ) and i ? n"' = Fn 6 '(w, b, a). 

The above generating function is a four variable generating function which is dependent 
on the size of the Aztec diamond. We could include a fifth variable marking the size of the 
Aztec diamond but that formula is less tractable. 

3.1. Boundary Generating function. Our method involves first finding a generating 
function for only those entries of K^ 1 whose row and column indices corresopnd to ver- 
tices on the boundary of the aztec diamond, called the Boundary generating function. This 
will be a generating function in three variables: two marking two location on the boundary, 
and a third marking the size of the aztec diamond. 

Consider an Aztec diamond with all edge weig hts equal to 1. Let Z n = 2 n ( n+1 )/ 2 which 
is the number of dimer coverings of an Aztec diamond of size n. For white vertices with 
Kasteleyn co-ordinates given by (2i + 1, 0) for < i < n — 1 and black vertices with Kasteleyn 
co-ordinates given by (0, 2j + 1) for < j < n — 1, let Z(i,j,l,ri) denote the number of 
dimer coverings of an Aztec diamond of size n, with edge weights given by 1, with the white 
vertex (2i + 1, 0) and the black vertex (0, 2j + 1) removed. 

Lemma 3.2. 

/„ -s Z(i,j,l,n) 1 ^ Z(i - k,j - 1,1, n - 1) 

^■") z =2 YZ^ %j)=(0,0),n>l 

n fc,/e{o,i} n 

(i_fcj_j)^(_i,_i) 

Proof. This proof refines the ideas in |27| which compute the edge placement probabilities 
for the Aztec diamond. 

Consider the Aztec diamond with the vertices (2i + 1,0) and (0, 2j + 1) removed for 
< i, j < re— 1. It is important to notice that the Kasteleyn orientation on the original Aztec 
diamond, restricted to this graph, remains a Kasteleyn orientation, so the corresponding 
entry of K^ 1 for the Aztec diamond can be computed from the partition function on this 
graph; since we have taken a = 1 for simplicity, this is just the number of matchings. 
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Figure 6. The dashed edges on the left are the two extra edges added as 
described in the proof of Lemma 3.2 for an Aztec diamond of size 4. The 
figure on the right is obtained from the figure on the left by applying urban 
renewal 16 times. 



The number of matchings is given by Z(i,j,l,n) and is in fact equal to the number of 
matchings of an Aztec diamond with two extra vertices vq and v± added with an edge between 
vo and (2i + 1, 0) and an edge between v\ and (0, 2j + 1) for < i, j < n — 1. We shall call 
vo and v\ auxiliary vertices and the corresponding edges auxiliary edges. To this graph, we 
apply urban renewal on the faces (2i + 1, 2j + 1) for < i, j < n — 1, i.e. we apply urban 
renewal n 2 times, see Figure [6] for an example. Each edge weight on the urban renewaled 
faces now reads 1/2 while there is an extra factor of 2 for each urban renewal step. We can 
remove all the edges connected to vertices incident to one edge due to them always appearing 
in each matching and apply edge contraction to the vertices incident to two edges. This gives 
an Aztec diamond of order n — 1 with two extra vertices Vi and v% where V2 has edges to the 
vertices (2i + 1,0) and (2i — 1,0) while V2 and v% has edges to the vertices (0,2j + 1) and 
(0, 2j — 1), provided that they are not corner vertices, see Figure [7J These extra edges can 
be seen from applying the removal of pendant edges to the auxiliary edges. In other words, 
we have 

(3.8) Z(i,j,l,n)= Z(i-k,j-l,l/2,n-l)2 n2 

k,ie{o,i} 

For the corner vertices under this urban renewal, we have Z(n — l,j, l/2,n — 1) = 0, 
Z(-l,j, l/2,n-l) = provided j ^ -1, Z(i, n-1, 1/2, n-1) = and Z(i, -1, 1/2, n-1) = 
provided i ^ —1 due to these shapes being not tileable. We have that Z(—l, —1, 1/2, n— 1) is 
equal to the number of dimer coverings of an Aztec diamond of order n—1 with edge weights 
1/2 and we have 2 n Z(—l, —1, 1/2, n — 1) = 2 n Z n _\ by performing the gauge transformation 
which multiplies all the edge weights incident to white vertices by 2. Therefore, we have 

(3.9) Z(i,j,l,n)= Z(i - k,j - l,l/2,n - l)2 n2 +I {lJ)={m Z n ^2 n . 

fc,ie{o,i} 
(i-fcj-O^(-i-i) 
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Figure 7. The result of edge contraction from Figure [6} The four dashed 
edges represent the result of the edge contraction on the two dashed edges in 
Figure [6j 




We can apply the gauge transformation which multiplies all the edge weights incident to 
white vertices by 2 to Z(i — k,j — 1, 1/2, n — 1). As there are n(n — 1) — 1 white vertices, we 
have 



(3.10) Z(i,j,l,n)= £ Z{i-k,j-l,l } n-l)2 n - l +l {ij)={m Z n ^T 

fc,/e{o,i} 
(i-k,j-i)^(-i-i) 

Dividing by Z n and noting that Z n = 2 n Z n _i gives the lemma. 



□ 



We can now compute the boundary generating function for removed vertices. Let 

oo n— 1 n—l 



(3.ii) zKM,*) = £££^f^ 



n=0 i=0 j=0 

Lemma 3.3. 

(1 - z)(2 — z(l + b)(l + w)j 



Proof. We can multiply (3.7) by z n and sum from < i < n — 1 and < j < n — 1 
which gives 

(3.13) Z(io,&,l,z) = -(l + &)(l+w)2Z(u;,M,z) + 



2 v -y 

Rearranging the above equation gives the result. 



□ 



We can now use the recurrence relations to find the generating function for the inverse 
Kasteleyn matrix. 



14 SUNIL CHHITA AND BENJAMIN YOUNG 

Proof of Theorem 3.1 for a = 1. The equations K.K~ l = I and K .K = I give two recur- 
rence relations. We have 
(3.14) 

K 1 (x+ei,y)5 Xl<2n +iK 1 (x+e 2 ,y)5 Xl>0 +±K 1 (x-e 2 ,y)5 Xl<2n +K 1 (x—e 1 , y)5 Xl>0 = 5 x=y 

where x,y G B, x = (xi,x 2 ),y = (2/1,2/2) and 

, , , J 1 if xi > 

(,3 - 15) ^ 1>0 " \ otherwise 

We also have for x, y G W 
(3.16) 

K 1 (x,2/+ei)^ 1 <2n+i-^ 1 ( a; 5 y+ e 2)<J J /2<2n+^ 1 (x, 2/-ei )<5j, 2 > +iif (x, 2/-e 2 )(5 ?/2 >o = ^=3,. 

We can multiply ( 3.14[) by w x b y for x, y G B and sum both quantities over B. Each term on 
the left-hand side of (3.14) can then be simplified to term involving G n (w,b) and a remainder 
term. For example, the first term of the left-hand side of (13.14) gives 



K 1 ( x + e i,y) S xi<2nW x b v 



1 



(3.17) 



K- 1 (x,y)w x b y 



G n (w,b)-£ i< xi < 2n K- 1 ((x l ^),y)w x l H 

(ai,0)gW,j/gB 

W\W 2 



Continuing the same procedure on the rest of the terms in (3.14) we obtain 

1 



(3.18) 



1 Wl w 2 , 

1 i H i + w\w 2 G n (w, b) 

W\W 2 w 2 w\ ' 



+ w ±i)G» n (w 1 ,b 1 M) 

W\W 2 w 2 



wiw 2 + — i ) 61, 62)^2 



2n 



•l -(wibi) 2n+2 \ / w 2 6 2 (l - (w 2 fr 2 ) 2n ) 
1 - W \b\ ) \ 1 - w\b\ 

where G\ (wx, 61, b 2 ) and G™(u>i, 61, b 2 ) are the coefficients of u>f> and u> 2n respectively in 
G„(w,b). 

It remains to compute G°(wi,6i,6 2 ) and GJJ(u>i, 61, 62)- Let 



(3.19) 



l#'(v,b) = 2 K-\{x x ,i),{j,y 2 ))w?wilP l bf 



l<xi<2n-l, 
l<y 2 <2n-l 



where i,jf G {0, 2ra}. We can compute G^(w\, b\, b 2 ) using the recurrence given in (3.16) and 
the method described above. We obtain 



4- + r 1 + 6162 + c°(u;i, 61, 6 2 ) - (j± 



+ gi]nr(-,b) 



i + 6i6 2 )^°' 2 >,b) 



1 - (wi&i) 1 



For G™(wi, 61, 6 2 ), using (3.16) we obtain 
6 



(3.21) 



^ + t) G ^ ~ {bk + 't) H ^ » 



i + hb 2 ) Hr>,b) = <6 2 ^ 1&l(1 " W " nb " n) 



1 - w 2 6 2 
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We can compute Hll 3 (w, b) using Lemma 
(3.22) H n {w,b) = 



3.3 



in the following way: let 



i 1 



(l+wl)(l+6l) 
2 



1 



(l+mi)(l+bi) 
2 



We have defined H n (w, b) to be the boundary generating function of the uniform Aztec dia- 
mond with the Kasteleyn orientation with the same co-ordinate labeling given to Z(w, b,l,z). 
This function is obtained by computing the z n coefficient of Z(w, b, 1, z) (a polynomial in w 
and b) and finding the appropriate sign by evaluating one of the terms. 
From H n (w, b) we can find H % n for i, j £ {0, 2n}. We find that 

(3.23) H°> 0f ~ ° - u <- 2 ^ 

(3.24) 



(w,b) = H n (w 1: 63)^162, 
H% 2n (v,b) = Hni-l/wl-b^w^b^hi, 



(3.25) 

and 

(3.26) 



< u (w,b) = H n (- W (,-l/bi) Wl wfb 



2ni 2n-l . 



#f' 2 >,b) = ^(l/^l/fc^f-^f 6f b 



2n-l„,.2n. 2rar 2n-l 



We can substitute (3.20) and (|3.21 ) back into ( 3.18| ) and use the expressions obtained 
from equations ( |3.23| > Q3.24D fl3.25j and ( |3.26p to find G„(w,b). 

□ 



Proof of Theorem 3.1 for all a (sketch). Generalizing this proof to arbitrary a requires no 
new ideas, but rather increases the complexity and amount of bookkeeping required. As 
such we have chosen to omit the details, but describe the main changes. 

For biased tilings, the boundary generating function and the recurrence relations from 
K.K' 1 = K~ l .K = I are modified. We also have that Z n = (1 + a 2 ) n ( n+1 )/ 2 . 

For the boundary recurrence relation, we can find the recurrence for the South West corner 
of the diamond, that is, removing the vertices (2i + 1, 0) and (0, 2j + 1) for < i, j < n — 1. 
We find that 

/Q o 7 a Z(i,j,a,n) ^ k+l Z(i-k,j -l,a,n-l) 

(3 ' 27) = ^ a Z n _!(l + a2) +I ^)=(°'°)- 

All boundary generating functions for K' 1 can be computed using manipulations of (3.27) 
and these are given in equations (3.1), (|3.2l), (3.3) and (3.4) 



The recurrence relations from K.K 
as described above. 



I and K .K = I can be used in the same way 

□ 



4. q vo1 WEIGHTING 

As mentioned in the introduction, it is possible to associate a discrete stepped surface to 
a tiling of the Aztec diamond. This function is called the height function. The most concrete 
way of viewing the height function is as the surface of a certain stack of blocks, called Levitov 
blocks |23l I24j . One can construct an edge weighting of the Aztec diamond in several ways, 
so that the contribution from each tiling is proportional to g#{ Levltov Blocks} _ adding 
a block to a tiling multiplies the weight of the tiling by q. Concretely, this means that for 
each face, the product of the edges parallel to e2 divided by the product of edges parallel to 
e\ is equal to q or g , depending on the parity of the face. Such a choice of weightings is 



16 



SUNIL CHHITA AND BENJAMIN YOUNG 



Figure 8. The two q vo1 weightings used in this paper for Aztec diamond of 
size 3. The figure on the left is q co1 weighting and the figure on the right is 
the r/ ,lim '. 




not unique, and we refer to any such choice as a g vo1 weighting of the Aztec diamond. For 
technical reasons, we make use of two such weightings here. 

We call the weighting q co1 to be the choice of weights where all edges parallel to e\ have 
weight 1 and the edges parallel to e<i have edge weights organized in columns given by 
aq 2n , aq l ~ 2n , aq 2n ~ 2 , . . . , aq -1 reading from left to right - see Figure [SJ We let K co \ denote 
the Kasteleyn matrix with q coi weights and its entries are given by 



(4.1) 



K co \(x,y) 



aq ~2n+x 1 -l i 

aq 2n ~ xi i 




y 
y 
y 



X 
X 
X 



±ei 
-e 2 



otherwise 



where x = (x%, x 2 ) € B and y G W. 

Theorem 4.1. The entries of the inverse of K co \ are given by 



(4.2) 



K-\(x,y) 



fl(x,y) ifx 1 <y 1 + l 

fi(x,y) + h{x,y) xi> + l 



for x = (xi,X2) £ W and y = (2/1,2/2) £ B, where 

j_(-xi+X2+yi-j/2+2)/2 JA+An-xi+X2 -y\-yi) (x\-X2 -j/i+?/2)/4 



h(x,y) 



(2tt±) 



(4.3) 



l/a,q 



dz dw 



w 



2/1/2 



n 



X2/2-1 

k=0 



( z + ag -2fc-l+w-l) lY^/ 2 -\azq 2k + x *+ 3 -yz _ x) 



Z (xi+1)/2( 1U _ z ) 
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and 



tt(*1-1/i-3)/2/ 2n+x 1 -4-2k\ r\(x2-y2-3) /2 . 2n-x 1 +4+2k 

f2(x,y) = ^ _ill^o ^ ; 



l/a,q 



1)2/1 / 2 -(, 1+ i)/ 2 nli / 2 - 1 (- + ag- 2 ^ 1 ^- 1 ) nn a/2 ~ W fc+a * +3 - w - 1) 



nl=o 1V V + aq 2k - 1 ) n^? 2+1)/2 («^ 2fc+2 - 1) 



where Tq is a contour surrounding the origin and ^y aq is a contour surrounding l/(aq 2 ), 
l/(aq 4 ),. . . , l/(aq 2n ) which does not intersect with Tq. 

When q = 1 , we recover the contour integral for domino tilings of the Aztec diamond with 
edge weights 1 and a for the horizontal and vertical tiles. 



4.1. Inverse Kasteleyn Entries on the North and East boundaries. We focus on 
computing the entries of the inverse Kasteleyn matrix for q co1 weights on the north boundary 
(the white vertices (2k — l,2n) for 1 < k < n) and the east boundary (the black vertices 
(2n, 2k — 1) for 1 < k < n). For this subsection, label the white vertices on the north 
boundary to n — 1 with the label i representing the vertex (2n — 2i — 1, 2n). Similarly, label 
the black vertices on the east boundary to n — 1 with the label j representing the vertex 
(2n,2n- 1 - j). 

Although we can find a boundary recurrence relation the q co[ weighting, it seems in- 
tractable to find an explicit solution. This inconvenience can be bypassed by working another 
(gauge equivalent) choice of q vo1 weighting. 

We let q dia s be the g vo1 weighting where all edges parallel to ei have weight 1 and all 
edges parallel to e<i on the North and East boundary have weight a which determines the 
rest of the edge weights - see Figure [8j Let if diag denote the Kasteleyn matrix for the g dia § 
weighting. Its entries are given by 



(4.5) K diag (x,y) = 

where x = (xi, X2) € B and y £ W. 

Lemma 4.2. The gauge transformation from q co1 weighting to g diag weighting is given by 
multiplying the weight of each edge incident to the white vertex i on the north boundary by 
q( t+1 } and the weight of each edge incident to the black vertex j on the east boundary by q 1 . 

Proof. Consider the g co1 weighting with the above multiplication of the edge weights. Fill 
in the remaining edge weights so that the horizontal edges have weight 1. The weighting 
recovered is exactly the q co1 weighting of the Aztec diamond. □ 

Let Zn iag (a, q) = Zn' ag (l; a, aq, aq~ 2 , . . . , aq 2n ~ l ) denote the partition function of tiling the 
Aztec diamond with g diag weights where 1 is the horizontal edge and a, aq, aq~ 2 , . . . , aq 271-1 
represent the vertical edges weights on the main diagonal (i.e. edges that intersect the line 
(0,0) to (2n, 2n)). Let Z diag (i,j,a,q,n) denote the partition function of tiling the Aztec 
diamond with g diag weights with vertices (2n — 2i — 1, 2n) and (2n, 2n — 1 — j) removed. 



1 y — x = ±ei 

± a q— 2mia(n— xi/2,n— (x2+l)/2) y _ x = e ^ 

iffl 2min(n-a;i/2-l,n-(a;2+l)/2))+l _ x _ 

otherwise 
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Lemma 4.3. 

(4.6) 



Z d ^(iJ,a,q,n) = 

7 diag - - ? 



Z n ia,g (a,q) 



J2 a fc+/ 2 ' diag ( i ~ - l,aq,Q,n- 1) + aI (i,j)=(o,o), 



n>l 



(fc,z)e{o,i} 

W(-i-i) 



(l + a^)Z n d !_ a f(ag,g) 



1 + a 2 <? 



/ 



Proof. We first find the recurrence relation for the partition function. By using urban re- 
newal, gauge transformation and removing the pendant edges, we find 



(4.7) Z^(a,q) = {l+a 2 q) n Z^ g {l-a,aq\aq- 2 ,aq\...,aq- 2n+ \aq 2n - 1 ) 

Therefore, we have 



(4.8) 



Zn &g (a,q) 



(1 + a 



2 n \n 



Z n-i(^ a( l ,aq ,...,aq ). 



at— >aq 



By ensuring that the local moves give the same change to the partition function and that 
the minimal configuration is the same, we find that 



(4.9) 



Z^(a,q) = (l+a 2 qTZ^f(aq,q). 



The rest of the proof follows from using the same procedures found in Lemma [3.2[ we add 
an auxiliary edge between an added auxiliary vertex and the vertex at (2n — 2i — 1, 2n) and 
another auxiliary edge between another added auxiliary vertex and the vertex at (2n, 2n — 
2j — 1). For this new graph, we can use urban renewal, gauge transformation and edge 
contraction. The factor a k+l given in the statement of the lemma is due to the reduction 
from i i — y i — 1 (and similarly j i— )■ j — 1) when n4n-l removes a vertical edge which 
has weight a. Similar to computing the partition function recurrence, we have to apply the 
weight change transformation given in (4.7) and (4.8) for the recurrence of Z dmg (i,j,a,q,n) 
but this changes the weight of the minimal configurations due to the removed vertices i — k 
and j — I for k, I G {0, 1}. To ensure that the minimal configurations remain the same under 
the reduction, we need to multiply by g t+J ' + . 

□ 



Let Z^ ol (a, q) denote the partition function of tiling the Aztec diamond with q co1 and 
Z col (i,j,a,q,n) denote the partition function of tiling the Aztec diamond with q co1 weights 
with vertices (2 n — 2i — 1, 2n) and (2n, 2n — 1 — j) removed. 

we can compute a three variable generating function for Z dmg (i, j, a, q,n)/Z 



From Lemma 



4.3 



and extract out the relevant coefficient for the size of the Aztec diamond. From this, we can 
extract out Z diag (i,j,a,q,n)/Zn ia,g using residue calculus. Using a rather non-elementary 
change of variables, we can find a double integral formula for Z diag (i, j, a, q, n) /Z diag . In the 
next lemma, we show that the double integral formula we obtained solves the recurrence and 
so we obtain a formula for Z col (i, j, a, q, n) / Z^° l . 
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Lemma 4.4. 

Z co \i,j,a,q,n) g (2+i+i) 2 -i 




dz dw 



Z-i (2vn) 2 

(4.10) 



{w - z)z n ~ i Y[lZt 1 {aq 2k - 1 + w) \^ k=0 {-l + <W fc+2 ) 
Proof. The gauge transformation from the proof of Lemma |4.2| gives 
/ 4 (m )2 +7 - Z diag (i, j,a,q,n) = Z col (i,j,a,q,n) 

V " / ? ^diag vcol 

Let 

f(i,j,a,q,n) = — / / dwcte 

(4 12) 



(w - z)z^ n^'V? 2 "- 1 + w) nLot- 1 + ^m 2k+2 ) 

We need to show that 
Claim 1. 

( 4 - 13 ) "dl^ = f(i,j,a,q,n) 



Proof of Claim. We show that f(i,j,a,q,n) satisfies the recurrence given in Lemma 4.3 
When i = j = 0, we have that the integral with respect to w has one singularity at l/(aq z ). 
We obtain 

f ( 00aqn) = ^f f n^Q 1 ^ 1 + g) ^ 



g 2 /• n»- J („„2*-J 
2?H 



r z™(a- i 9- 2 -^)re=o 1 (i+«v fc+i ) 



(4.14) 



2?™ 



F z n (a- 1 ^ 2 - z) rifelo l 1 + « 2 9 2fc+1 



1 ' mZliaq 2 "- 1 + z) aq 



2ni J Vo z n-i (a-ig- 2 - z ) rK=o(l + a?q 2k + l ) 1 + a?q 

where the fourth equation follows from pushing the integral with respect to z through infinity. 
By taking the change of variables w h-> q~ l w and z i-> we can compute /(0, 0, ag, q,n — 
1). We obtain 

(4.15) 

w n - l \[ n k =\(aq 2k - l + z 



f (0, 0, ao, q, n — 1) = rr / / = dw dz 

V ' (2vrz) 2 7r A-, ^"H™ " z) ^(ag 2 *- 1 + w)(waq^ - 1) 



i f nLo( a ^ 



z 



2ni J To z n-l (a-lg-2 - *) n£j(l + a 2 ^ 1 ) 
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We can compare equations (4.14) and ( |4.15| ). We obtain 
(4.16) /(0,0,a,g,n) 



1 + a 2 q 1 + a 2 

When i,j ^ 0, we can write 

f(i,j,aq,Q,n) = J^) 2 / / dwdZ 

(4.17) r ° T ^ q 

w n ir k Zoaq- 2k+2{n - j - 1) + z 



z n _ z ^ Y[^ = q 1 aq 2k + w nl=o a Q 2k+3 w — 1 
By manipulating the integrand of the above expression, we can write 

(2+2i+j)(l+j) r r 

f(i,j,aq,Q,n- 1) + f(i- l,j,aq,q,n- 1) = ( 2 7d) 2 ' ' 



To JT 

n-l Tin-2 



n-l 

aq,q 



z 



llfc= ag , 2(1+/) --1 



(4.18) = Xlft=u P — l + a^™z 



w — z) rife=o a{ ? 2fc + w ni=o a Q 2k+3w — i 



(2+2i+i)(l+i) /■ r ^-l jjn -1 2 fc+2(n-i-2) + J? 



(2vri) 2 7 ro 7 r n-i z n-i( w _ z ) 2 aq 2k + w -qj =o ag 2fc+3 w _ i 

We can also write 
(4.19) 

7 (2+2i+j)(l+j) 
(27Ti) 



f(i,j,aq,q,n- 1) + f(i,j- l,aq,q,n- 1) = 77 — g / / <iu;cfe 




gii—i—l i 



L aq,q 



u>»- i na°g- 2fc+2(w - j - l) +^ L- 2 (i+^) + «t 2(w+2 ^ 

_ g (2+2i+j)(l+j)-2i-2j-2 r r 112=0 ag- 2 ^ 2 ^'- 1 ) + Z 

(2vri) 2 " ir hill ^{w - z) l^Zf 1 a Q 2k + w Wi=o aq 2k+3 w - 1 

For the above equation, we can make the change of variables w h+ g 2 w and z i-+ <7 2 z. We 
obtain 

aq (2+2i+j)(l+j) 

J {'hJi^lili 'i' - -U t- J V^J ~ t,uq,q, n - ±| = 

(4.20) 



f(i,j,aq,q,n- 1) + - l,aq,q,n- 1) = / / . 



/ -pi / Tin — 1 
' 1 J 1 aq,q 



w"" 1 UkZl gq- 2k + 2 ^-i- 2 ) + z 
z n ~ l (w — z) rifc=-7 2 a 9 2k + w nl=i aq 2k+3 w — 1 
We can add ( |418| ) and ( [420] ). We obtain 
(4.21) 

f(hj,aq,q,n- 1) + /(i - l,j,aq,q,n- 1) + f(i,j,aq,q,n- 1) + /(i, j - l,ag,?,n- 1) 
ag (2 + 2 i+ , )( i + /) , , ^ (n^o 1 ag- 2 *+ 2 (^- 2 ) + z) (l + ^ffi ) 



(2Tri) 2 JtJt2-i z n -\ w - z) rifcj 2 v fc + ^ ni= a ^ 2fe+3 ^ - 1 

ag (2 + 2 4+J )(i + ,) (1 + a 2 {?) , , «,« (n*=5 ag- 2fc + 2 ^'- 2 ) + z) 



(2vri) 2 7 ro 7 r n-i ^n-i( w _ z ) nL-i 2 a « 2fc + w ni=o V fc+3 w - 1 



dw dz 
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Finally, we can take the change of variables w h-> q 1 w and z i— > q z and multiply by g i+J " +1 
which gives f(i,j,a,q,n) as required. 

□ 



Lemma 4.4 now follows directly. 

□ 



Proof of Theorem J^.t. Using Lemma 4.4 we can write K~ l for the entries on the NW bound- 
ary using the Kasteleyn co-ordinates and adding the relevant sign from the Kasteleyn orien- 
tation. We obtain 

(4.22) 

i (4n-a;i-y2+2)/2^(4+4n-2:i-?/2)(^i+J/2-4n)/4 r r 

K~ ((xi,2ri), (2n,y 2 )) = , ,, 2 / / dzdw 

z (*l+l)/2( u; _ z ) n ^-l)/2 (u , + aq 2k-l) n^ 1+1)/2 (mV fc + 2 - 1) ' 

From the above expression, we first move the white vertices to x G W by using the recur- 
rence generated from K.K -1 = I when no white and black vertices are adjacent under the 
recurrence. That is, we use the equation 

K- l {x,y) = -K- X (x + 2e 1 ,y)S Xl<2n -i- a±q Xi ~ 2n K~ x {x + e x + e 2 ,y) 

-aiq 2n - Xl - 1 K- 1 (x + e 1 -e 2 ,y)5 Xl<2 n-i 

for x = (xi,x 2 ) G W with x\ > and y G B. We can then move the black vertices to y G B 
using the recurrence K~ x .K = I provided no white and black vertices are adjacent under 
the recurrence. That is, we can use the equation 

K- l (x,y) = -K- 1 ( y x,y + 2e 1 )5y 2<2n ~i - a±q 2n ~ Vl K~ 1 (x,y + ei + e 2 )5 y2<2 n-i 
' Ul) -ai^ 1 - 2n+1 K- 1 (x,y + e 1 -e 2 ) 

for y = (2/1,2/2) G B with y 2 > and x G W. These two recurrences give fi(x,y) defined 
in ( |4.3[ ). If the vertices cross under the recurrence, then there is an extra factor counted 
from the contribution of z = w in the expression of fi(x,y). This is exactly the expression 
h{x,y). 

□ 



5. Two-Periodic Weighting 

In this section, we compute the generating function of the inverse Kasteleyn matrix for the 
two-periodic weighting of an Aztec diamond of size Am. This is much akin to the parameter 
a introduced in Section [3j except now there are two parameters, a and b, which decorate the 
edges of the aztec diamond in a checkerboard fashion described below. As remarked in the 
introduction, one of the special cases of this model is equivalent to a different, uniform tiling 
problem: namely, the so-called diabolo tilings of the fortress introduced in [27]. The dual 
graph of the fortress is the square- octagon lattice with certain boundary conditions. This 
model was the main motivation for this work. 

Because of this periodicity, the problem becomes complicated in two ways: the recur- 
rence relation increases in order, and the generating function becomes a vector of generating 
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Figure 9. The two-periodic weighting of the Aztec diamond with the pic- 
ture on the left being an even ordered Aztec diamond and the picture on the 
right being an odd ordered Aztec diamond. The edges surrounding each face 
are given the same weight, denoted by the parameter in the center of each 
face with the exception of the boundary of the odd ordered Aztec diamond, 
where the edge weights are given by the adjacent parameter. 




functions. Finally, the relation which we solve in order to compute the boundary generat- 
ing function involves a certain matrix multiplication which must be explicitly diagonalized 
before the recurrence can be solved. 

As the expressions are somewhat complicated, we can only write the overall generating 
function as linear combinations of the entries of the if _1 (w, b) where w and b are white and 
black vertices on the boundary of the Aztec diamond where the matrix K is given below. 
These entries of the inverse Kasteleyn matrix are computed explicitly and, for the bottom 



right hand boundary, are given in Lemma 5.1 



We consider a two-periodic weighting of the Aztec diamond and compute its boundary 
generating function. Due to the periodicity or the weights, we have two types of white and 
black vertices. We denote for i G {0, 1} 



(5.1) 

and 

(5.2) 



{(xi, X2) : (xi + X2) mod 4 = 2i + 1, (x\, x?) G W} 



Bj = {{xi,x 2 ) : (xi + x 2 ) mod 4 = 2i + 1, (xi,x 2 ) € B}. 



For an Aztec diamond of size n, we give weights a and b to the Aztec diamond in the 
following way: if the size of the Aztec diamond is even, i.e. n = 2r, then the edge weights 
around each face (2i + 1, 2j + 1) for 0<z,j<n — 1 are given weight a if (i + j) mod 2 = 
and weight b if (i + j) mod 2 = 1. Conversely, if n = 2r — 1, then the edge weights are 
obtained from embedding the diamond in an Aztec diamond of order 2r. Figure [9] shows 
this choice of edge weights. 
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Let K denote the Kasteleyn matrix for an Aztec diamond when n is even and K denote 
the Aztec diamond when n is odd. We have 



(5.3) 



and 



(5.4) 



K(x,y) 



a(l — i) + bi if y = x + e\, x G Bj 

(ai + 6(1 — i))i if y = x + e2, x £ Bj 

ai + 6(1 — i) if y = x — ei, x G Bj 

(a(l — i) + 6i)i if y = x — e2, x G Bj 



K(x,y) 



ai + 6(1 



if y = x + ei, x G Bj 



(ai + 6(1 — i))i if y = x + e2, x G Bj 



a(l — i) + bi if y 
(a(l — i) + 6i)i if y 



x 
x 



ei, x G Bj 
e 2 ,x G Bj 



for x G B and y G W. 

We now give the entries of for white and black vertices on the bottom and left 
boundaries respectively for K is defined in (5.3) and the size of the Aztec diamond is even. 
To shorten the length of the formulas, we will write [{[2 = i mod 2. 

Lemma 5.1. For an Aztec diamond of size n, let n = 4m and let K denote the Kasteleyn 
matrix given in (5.3). Let L(a,b,i,j) denote K ((2i + 1,0), (0,2j + 1)) for the Kasteleyn 
matrix with parameters a and 6. We have 



L(a,b,i,j) 



(5.5) 



(2vri) 2 

m— 1 

£ E 

r=0 fc,ie{0,l} 



\z\=U\w\=l 



U,Li/2j+l z L772j+l 



where for 1 < i, j < 4 we have gj(a,b,w, z) = Njj(a, 6, if, z) which 
and 

(5.6) a k (a,b,w^ - 



dw dz 

is given in Appendix \a\ 



for k G {0, 1} with 
(5.7) Pi(a,b,w) 



Arabia 2 + 6 2 )(^ v /(6 4 + a 4 )w + a 2 6 2 (l + w 2 )) fc 
(a 2 + 6 2 )^ - (-1) V(^ 4 + a 4 )^ + a 262 ! 1 + ™ 2 



^2a6(a 2 + 6 2 ) 



Note that the expressions Qo(a, 6, u>) and ai(a, 6, w) are polynomials in to and the L(a, 6, i, j) 
is a rational function in a and 6. 

The boundary inverse Kasteleyn entries have periodicity four with the size of the Aztec 
diamond. It is also possible to find the entries of the boundary inverse Kasteleyn matrix for 
the size of the Aztec diamond equal to 4m — 1, 4m — 2 and 4m — 3 but we will not state 
these results here. The proof of the lemma is given in the next section. For the generating 
function of K" 1 , we need the following terms: 

cgf(wi, W2) = 2(1 + a 2 ) + a(w\ + w^ 2 )(w\ + w^ 2 ), 

/ \ f 9 .9 9 9\ 9 9 9' 

^i,o{wi, W2) = — a (lOj u> 2 + w 1 w 2 ) — aiw-t + aiwi —2a 



(5.8) Sifi[wi,W2) = —a(w 1 u> 2 +w 1 u> 2 J - aiw 1 + aiw 1 

Si,2n{wi,W2) = w 2n ^— a {w\w2 + W\ 2 ^|) + o,i.w\ — ai.w^ 2 — 2a 2 ^ 1 ~ i ^ . 
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the sum of a geometric series and we also 



for i g {0, 1}. We set f n (x) = (1 - x n )/(l - x 
let for n = 2r and w = (wi, w 2 ) and b = (61, 62), 

fr{w A 2 b A 2 )f r {w\b\)b 2 wi{a + w\b\ + w\b\{\ + aw\b\)){w\ + b\- i(l + w^)) 



(5.9) 



d(w, b) 

+ ((1 - iwl) Wi%(l + aw 2 2 bl) + {wl - i) wl n+1 bl n b 2 (a + w 2 2 b 2 2 )) f r (w 4 2 b 4 " 

-Mwtbi 



^sides( w ) b) f / ,_ 4 4. 



CGF (&1,&2; 



where 



(5.10) 



4ides(w,b) = s 0fi (wi,w 2 ) [a(b\ - i)w x b 2 + w±b 2 1 (l - ibj)) 

+ s lfi (w 1 ,w 2 ) {{b\ - i)wfb 2 1 b 2 + atof&ifc^l - i&?)) 

+ s ,2n(wi, ^2) ((1 - i^Jwi^ 1 + a Wl bl n+ \l - xb\)) 

+ si,2n(^l, w 2 ) (o(l - i^K^"" 1 + ^6^+1(1 _ i6 2)) . 



We denote the generating function of K 1 for if defined in (5.3) as 
(5.11) G(a,6,w,b) = ^^if- 1 (x,y)w^ 

sGW J/6B 

where w x = w^w^ 2 and b y = b y ^b y 2 for x = (xi,x 2 ) and y = (3/1 , ^2) ■ 

Theorem 5.2. For an Aztec diamond of size n, the generating function of K~ x for K 
defined in (5.3) with b = 1, n = Am for m £ N, w = (ioi, w 2 ) and b = (61, b 2 ) is given by 

(5.12) 

d(w, b) 



G(a, 1, w, b) 



+ 



where 



cgf{w\-,w 2 ) 

v E 



E ^ ^ ^ 

vijS{0,l} fc,Ze{0,2n} (xi,fc)eWi (l,y 2 )&Bj 



(5.13) 



K- 1 ((x 1 ,0),(0,y 2 )) 
K- 1 ((x 1 ,0),(2n,y 2 )) 

^((x^n),^)) 
if- 1 ((xi,2n),(2n,y 2 )) 



:,k(wi,W2)sj,i(b 2 ,bi)K 1 ((x 1 , k), (£,3/2)) 
c GF (wi,w 2 )c GF (oi,o 2 ) 

xi - 1 2/2-1" 



L I a, 1 



2 ' 2 
i2n -l-x 1+y2L ^ a n . 

i 2n-l+*x-« 2Z/ 

Si + 1 



Xl + 1 2/2 ~ 1 

2 ' 2 

2/2 + 1 



L I a, 1, n 



-,n 



2/2 + 1 



and 



and L(a,b,i,j) is given in Lemma 5.1 



For th e pro of of this theorem, we are not able to follow the approach exactly as given in 
Theorem 3.1 because the relations K.K^ 1 = I and K~ 1 .K = I, while technically sufficient, 
are not of a suitably nice form to allow any progress. Instead, we need to use two further 
recurrences K* .K.K~ l = K*l and K~ l .K.K* = l.K* where K* denotes complex conju- 
gate transpose. These identities have an interpretation in terms of the discrete Laplacian 
interpretation, see [20] — though this interpretation is only heuristically relevant here. 
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5.1. Boundary Generating Function. In this section, we find the boundary recurrence 
relation and solve the recurrence. The method is similar to the one used in Section [3] but 
due to the two-periodic weighting, we have a different recurrence for each type of white and 
black vertex. This means we have a matrix equation explaining the recurrence which is also 
periodic. Due to the nature of the computations, we had to rely heavily on computer algebra 
in this subsection. 

Let Zp(a,b,n) denote the partition function of an Aztec diamond of size 2n with edge 
weights a and b described above and let Zp(a, b, n) denote the partition function of an Aztec 
diamond of size 2n — 1 with edge weights a and b as described above for an Aztec diamond of 
size 2n — 1. Let Zi-i(i,j, a, b, n) count the number of weighted tilings of an Aztec diamond of 
size 2n with vertices (4i + 2k + 1, 0) and (0, 4j + 21 + 1) removed with edge weights a and b. 
Let Zki(i,j, a, b, n) count the number of weighted tilings of an Aztec diamond of size 2n — 1 
with vertices (4i + 2k + 1, 0) and (0, 4j + 21 + 1) removed with weights a and b. 

We introduce the following generating functions: denote 

OO TL— 1 U— 1 ry I ■ ■ 7 \ 

Z k i(i,j,a,b,n) 



(5.14) 

and 
(5.15) 

for k,l 6 {0,1}. Let 



G H {a,b,x,y,z) 



n=0 i=0 j=0 



Zp(a, b, n) 



aula, 6. *, = %# »W. 

n^^i Z P {a,b,n) 

n=0 «=0 j=0 ' 



(5.16) G(a,6,x,y,z) 
and similarly, let G denote the corresponding vector for Gki- 



( G 00 (a,b,x,y,z) \ 
G i(a,b,x,y,z) 
G 10 (a,b,x,y,z) 
\ G u (a,b,x,y,z) J 



Lemma 5.3. The boundary generating functions G and G satisfy the following recurrences 
(5.17) 
and 

G(c, d, x, y, z) = zC(c, d, x, y).G(c, d, x, y, z) + D(c, d)- 



G(a,b,x,y,z) = A(a,b,x,y).G ( — , — , x, y, z ) +B(a)-^— 

2a 2b j 1 — z 



(5.18) 
where 

(5.19) 



A(a,b,x,y) 



1 


2/ 


a; 


4j 2 




4a 2 

X 




lab 

y 


4ab 


4^5 




AaT) 


iP 


iP 


W 



\ 



B(a) 



1 - z 






V o J 



(5.20) C(c,d,x,y) 



+ d 2 



/ d 2 

d 2 
d 2 

V 



cdy 
cd 
cdy 
cd 



cdx 
cdx 
cd 
cd 



c 2 xy 

2 

C X 

c 2 y 



and D(c,d) 



c 2 +d 2 
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Proof. By urban renewal, edge contraction and removal of pendant edges, we can find that 
Zp(a,b,n) = (4ab) 4n Zp(l/(2o), 1/(26), n). We now write a relation for Z st (i, j, a,b,n) for 



s,t £ {0,1} following the same method given in Lemma 3.2 That is, we add two auxiliary 
vertices to the graph x>\ and v<i and their two auxiliary edges which are edges from v\ to 
the vertex at (4i + 2s + 1,0) and from V2 to (0,4j + 2t + 1) so that the resulting graph 
is planar. We can then apply urban renewal, edge contraction, removal of pendant edges. 
When splitting the additional square from the auxiliary edge obtained from v±, we have an 
extra factor of l/(2a 1_s 6 s ) for s £ {0, 1} obtained from the edge weights of the additional 
square. We have a similar factor for the splitting the additional square from the auxiliary 
edge corresponding to V2- We find that for s,t £ {0, 1} 

(4a6) 4 ™ 2 / 1 1 \ 

(5.21) Z st (i,j,a,b,n)= -^—^Zm (i - k(l - s), j - 1(1 - t), — , -, n J . 

fc,«e{o,i} 



The boundary conditions for Zj~i are similar to those given in Lemma 3.2 We find that 



(5.22) 

Z st (i,j,a,b,n) 



E 

(4a6) 4n2 



(4a6) 



An z 



(i-fej-O^(-i-i) 



4a 2 - s -*6 s +*' 



Z« \ i-k(l-s),j-l(l-t), 



1 1 

2a' 26' 



, n 



2a 



Z P 



1 1 



2a 2b 



'^ n ) I (i,j>,t) = (0,0,0,0) 



We can divide both sides of the above equation by Zp(a, 6, n) which gives 
(5.23) 



Z st (i,j,a,b,n) 
Zp(a, 6, n) 



E 



Zfcz (i - fe(l - s),j - 1(1 - t), i,n) I(ij>,t)=(o,o,o,o) 



fc,;e{o,i} 



tea— *6^Z P (i i,n) 



+ 



2a 



For the recurrence equations given in the equation above, we can multiply by x % y d z n and 
sum over the relevant quantities. This gives a system of equations involving the generating 
functions whose matrix equation is given by (5.17). 

Analogously, we can find a recurrence relation for Z s t(i,j,c,d,n). We find that for s,t £ 
{0,1}, 



(5.24) 

Z st (i,j, c, d, n) 
Zp(a, 6, n) 



E 



•2-k-l„k+l 



c k+l Z kl (i - k(l - s),j - 1(1 - t), c,d,n- 1) 



fe,/e{o,i} 
(i-kj-i)^(-i-i) 

c 



c 2 + d 2 



Zp(c, d, n — 1) 



+ c 2 + d2 Hi,j,s,t)=(0,0,0,0)- 



The above recurrence relation gives the generating function given in (5.18). 

Let 
(5.25) 

M(a, 6, x, y) = A(a, 6, x, y).C ((2a)- 1 , (2b)-\x, y) A {(2a)-\ (2b)-\x, y) .C(o, 6, x, y) 



□ 
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and define the vectors 

(5.26) Bi(a, b, x, y) = A(a, b, x, y).D ((2a)- 1 , (26)" 1 ) + B(o) 
and 

(5.27) B a (o, b, x, y) = A(a, 6, x, y).C ((2a)- 1 , (2b)-\x, y) .Bx ((2a)" 1 , (2b)~\ x, y) . 

Lemma 5.4. For n = 4m, the n th coefficient of z of the generating function G(a, b, x, y, z) 
is given by 

(5.28) 



m— 1 



^rAT- 1 .(B 1 + B 2 ) 



i=0 

where (A,T) is the eigensystem of M(a,b,x,y). Explicitly, the eigenvalues of M(o, 6, x,y) 
are given by 

(5.29) X 2 i+j+i(a, b,x,y) = Pi(a,b, xf f3j(a,b,y) 2 

and the eigenvectors are given by 
(5.30) 

a?b 2 (l +x)(l+y) ' 
((P-a^x-i-iy^Vi^+b^x+^^il+x 2 )) 
ab( l+x) 

(b 2 -a 2 )y-(-iy^/y^(ai+bi)y+a 2 b 2 (l+y 2 ) 
ab(l+y) 
1 



V2i+j+l 



\ 



) 



Although the above lemma does give the boundary generating function for a two-periodic 
Aztec diamond, the expression is complicated. We believe that the expression given in 
Lemma 5.1 is more feasible for potential asymptotic computations. 



Proof. From Lemma 5.3, we have a generating function equation for G and G. We can write 
G(a, b, x, y, z) = zA(a, b, x, y).C ((2o)" 1 , (2b)-\x, y) .G ((2a)" 1 , (26)" 1 , x, y, z) 

'"" Ul + (A(a,b,x,y).B ((2a)-\(2b)- 1 )+B(a)) 

Apply the above equation to itself, we can write 
(5.32) 

G(o, b, x, y, z) = z 2 M(a, b, x, y)G(a, b, x, y, z) + Bi(o, b, x, y) + — -B 2 (a, b, x, y) 

= z 2k (M(a, b, x, y)) k . ( ^^1^, b, x, y) + ^B 2 (a, 6, x, y) ) 

k=o \ z z J 

by using the expansion of a geometric series of matrices. The above equation can be solved 
but it does not appear to give a tractable answer. As we are interested in the n th coefficient, 
we can use the above expansion and the expansion of (1 — -z) -1 to find the coefficient of z n . 
This is given by 

m— 1 

(5.33) ]T (M(o, b, x, y)) k . (B x (a, b, x, y) + B 2 (a, b, x, y)) . 

k=0 

The eigenvalues and eigenvectors of M can be verified by M.u = Xv where A is the eigenvalue 
for the eigenvector v. □ 
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We can now prove Lemma 5.1 



Proof of Lemma 5.1 . From Lemma |5.4[ we can write the n th coefficient of G(a, b, x, y, z) as 

m— 1 4 

(5.34) EE^j 

where {Xi}f =1 are the eigenvalues of M(a, b, x, y) and Xj are four column vectors which are 
the coefficients of in the following expression 

(5.35) T.diag (ax, a 2 , a 3 , 04) T 1 . (Bi(a, b, x, y) + B 2 (a, 6, ac, y)) . 

where diag(ai, 02, 03, 04) denotes a diagonal matrix with four entries 01,..., 04. We can 
rewrite (5.34) as 

m— 1 

(5.36) 



A;=0 



where 



(5.37) 



/ 1 


1 


1 






f 


K \ 




1 - 


1 


1 


-1 






4 




1 


1 - 


1 


-1 






a 3 




V 1 - 


1 - 


1 


1 J 






A 4 J 










( 1 


1 


1 


1 \ 


2 X3 


x 4 




1 


-1 


1 


-1 


]■ 


1 


1 


-1 


-1 








V 1 


-1 


-1 


1 / 



A computation shows that for 
(5.38) 

/ 



Di = 16a6(a^ + 6 2 ) 2 diag 



that 



1 



Vj/ v /(b 4 + a 4 )y + a262(i + y 2) 

^\/( fe4 + a ) x + « 2ft2 (! + 
\ ^\A ft4 + a 4 )^ + a 2 ^ 2 (l + x 2 )Vy^/{b A + a 4 )y + a 2 ^ 2 (l + V 2 ) J 



N(a,b,x,y) = Y.Di 



(5.39) 

where N(a, 6, x, y) is the matrix defined in ( A.l ) and 

( al{a,b,x)al(a,b,y) \ 
a l (a,b,x)a\(a,b,y) 
a\(a,b,x)a l (a,b,y) 
\ a\(a,b,x)a\(a,b,y) J 



(5.40) 



Di 



V 



1 \ 



/ 



V A| / 
This means that the 



where ay is defined in the statement of Lemma 5.1 for j 6 {0, 1}. 
coefficient of z n in the expression G(a, 6, x, y, z) is equal to 

/ a (a,b,x)a l (a,b,y) \ 

(5 41) N(abxv) ^(a,b,x)a\(a,b,y) 

^ > r*{a,o,x,y). ^ a ^ x )^ a ^ y ) 

\ a\(a,b,x)a\(a,b,y) J 

We can extract the relevant coefficient of the above equation to find \L(a,b,i, j)\. As each 
matching of the Aztec diamond with the vertices (2i + 1, 0) and (0, 2j + 1) removed has the 
same sign, the sign of one such matching determines the sign of L{a, b, 
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□ 



5.2. Generating Function of K . We now prove Theorem 5.2 

For notational purposes in the proof, we write G = G(a, 1, w, b) and let 



(5.42) 



G\ 



x=(i,j) 



K- 1 (x,y)w x by 



X 1 =l,X 2 =J 

J/SB 



for i G {1, 2n — 1} or j £ {0, 2n}. For j £ {0, 2n}, we also write 



(5-43) G\ x=(xuj) 

and for i 6 {1, 2n — 1} 

(5-44) G\ X=M - 



E G| 



a;=(i,i) 



l<i<2n-l, 
i mod 2=1 



E G \ 

0<j<2n, 
i mod 2=0 



x=(i,j) ■ 



Proof of Theorem 5.2. Recall that K* is the conjugate transpose of K. Applying K* to both 
sides of the equation K.K^ 1 = I, we obtain 

(5.45) 

a(K^ 1 (x + 2ei,y)5 X2<2n o~x 1 <2n-i + K^ 1 (x + 2e 2 , y)^x 2 <2n^x 1 >i 
+ K' 1 (x - 2ei,y)5 X2>0 5 Xl> i + K~ 1 (x - 2e 2 ,y)5 X2>0 5 Xl<2n -i) 



+ 2 1 + a 2 



J X 2 : 



=o(^xew + ^^x&xj 



Jx2=2n 



[a 5 x£Vo + 5 xeVl )) K (x,y) 



+ a±S X2=0 (-K 1 (x + e 2 -ei,y)5 Xl>1 + K 1 (x + ei - e 2 , y)^ 1 < 2 n-i) 

+ ai5 X2=2n (^ _1 (x + e 2 - ei,y)5 Xl> i - K~ 1 (x + ei - e 2 , y)(5 Xl<2n _i) = i^*(^ = ^(x)) 

for x = (xi,x 2 ) G W and y G B and where i£"*(<5s =2/ (:c)) means JT* is applied to the white 
vertex y when it equal to the vertex x. We also have a similar expression for a recurrence 
that moves the black vertices. 



We can multiply both sides of equation (5.45) by w x h v = w^ 1 2 fef 1 6 2 2 and sum over all 
the white and black vertices of the Aztec diamond. We can then rewrite each term of (5.45) 



using G(a, 1, w, b), expression with white vertices on the boundary and expression with x\ = 1 
or ii = 2n — 1. For example, we have 



^Vb^if 1 (x + 2e 1 ,y)5 X2<2n 5x 1 <2n- 



(5.46) 



G(a, l,w,b) - ( Y,xi=l,x 2 ^o + Yl £2=o 



2„„2 



W? W 



l w 2 
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We can simplify the rest of the terms of (5.45) and combine them. We obtain 
(5.47) 

(2(1 + a 2 ) + (wj + w^ 2 ){w 2 2 + w 2 2 ))G 
-aG\ ( 1 + wA A-aG\ + < 

1 + Wo Wn , \ _ ( 1 + Wn 1 



+ a G \x=(lfi) ^^2" -^2- W l i ) +a G U(l,2n) ^^2" " ^2 + »1 1 

^yi /^(l + Wo) 2 2 2 \ ^-yi /W(l + ^2) W ? 

+ a G| x=(2n _ 10) - u, 2 ^ 2 + t«?i ) + a G| x=(2 „_ li2n) ( ^ - J - ™ 2 i 



2. 

... _ .. , , ; - 

W% J y ' V Wo Wo 

2„..-2 , „..2„..-2\ „„„2- 1 „„,.-2j 



+ G| x= (a; 1( Q) (— a(u; 1 u> 2 + w\w 2 ) — aw l ± + aw 1 i — 2) 
+ G| a . = ( a . lt o) (— a(u>^ 2 u>2~ 2 + u^u^ 2 ) — aw 2 i + awj" 2 i — 2a 2 ) 

a;eWi 

+ G\x=(x u 2n) {-a(wjwj + w- 2 u>|) + aw\l - aw^ 2 ± - 2a 2 ) 

rr6W 

+ ^U^a^n) (-«(^i^2 + w^wj) + awii - aw^ 2 i - 2) = ^ K*(4 =J/ (x)) 

The terms G| 2 . = ^ lx2 ^ and G\ x= , 2n _ lx \ i nv °l ve white vertices away from the boundary. 
We can use the recurrence relation from K.K -1 = I to write these expressions in terms of 
the boundary vertices. For G\ x= ^ X2 y we use the recurrence relation 

(5.48) K^ 1 (x + e 2 + ei, y) = -^—5 x+e2=y - iK'^x, y) 

for x = (1,2:2) S Wi for i £ {0, 1}. We can multiply the above equation by w x b y and sum 
over the black and white vertices. We find that (5.48) becomes 
(5.49) 

w 2 2 G \x=[i,x 2 )~ w l G L=(i,o) = Yl 2 -iZiS x+e2=y w x b y -± G| x=( i )0 ) + i G\ x 



«e{o,i} xi=i 

X2f^2n,x£Wi 



,-(l,2n)- 



We can write 



(5.50) J2 E ^x + e 2 = y w x by = f n/2 (wibi)wMa~ l + 



, _ . . 

a x 

«e{o,i} x i =1 

X2^2n^x&ti 



by evaluating each sum separately. Using ( 5.49| ) and the above equation, we find that 
(5.51) G\ X={1:X2) (w 2 2 + i) = L/2{™ibt)wMa- 1 + w 2 b\) + w 2 2 G\ x={1>0) + i G\ x=(1>2n) 
We can do a similar procedure to find G f | x _/ 2n _ 1 x s using the recurrence 



(5.52) K 1 (x + e 2 + ei,y)i = -jd x+ei=y - K 1 {x,y) 
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for x = (2n — 1,2:2) G Wj for 2 £ {0, 1}. By performing a similar procedure as above, we 

obtain 

(5.53) 

G| x =(2n-1,* 2 ) K~ 2 i+1) = /n/2(^ 4 ^X- 1 6? n 62(l+^ 2 ^a- 1 )+t, 2 - 2 i G^n-l^U^-l 

We can substitute equation ( 5.51| ) and (5.53) back into (5.47). We find that the coefficient 
of a G f | a!= fi 0) m (5-47) is given by 

(5.54) 1 + ^ w 2 



wj(l + w 



W7 



[w 2 2 + i)w 2 w% 



0. 



Similarly, we also find that the coefficients of a G\ x= ^ 2n \, a G\ x , 2n _ 1 ^ and a G\ x= ^ 2n _ 1 2n > 
are zero. This means we can reduce ( |5.47 ) to 



c G f(wi,w 2 )G 



a(l + w 2 )w\b 2 (a 1 + w 2 b 2 



fn/2(btwp 



wlwl(w 2 



+ i) 



(5.55) 



a(l + w|)u; 1 n+1 6f6 2 (l + w^a" 1 ) 



Wo\W n 



+ 



fn/2(b 2 W 2 ) 



^2 s ifi (wi,w 2 ) G\ x={xu0) + 
ie{o,i} x ' eWl 



?:e{o,i} 



Si,2n(wi,W 2 ) G\ x= ( xi:2n ) = ^ K* (5( x=y) (x)) . 



To evaluate iT* (5 x=y (x)) we have that 
(5.56) 



K*(8 x=y (x)) 



Jx 2 >0 



a* 5 : 



x-ei=y 



ia % 5. 



x-e2=y 



)+S.. 



X2<2n 



- , x+ei=y 



10, 



We can multiply both sides by w x b y and we obtain 
(5.57) 

^wWT(^ =y (x)) = 



5^ ^(ttJibi)* 1 ^^; 

ie{o,i} zew; 



•A'2 



a' 

0.7) 



-a 1 



.''2 



>0 + 



61b: 



-£x 2 >0 + a % b\b2&x2<2n 



l-i- b2 x 
a 1 ~[~ X2<2n 
01 



This can be evaluated explicitly and gives the first expression of <i(w, b). Hence we obtain 
an expression for G(a, 1, u>, b) as linear combinations of K ((x%, 0), y) and ivr _1 ((xi, 2n), y)) 
for 1 < xi < 2n — 1 and y £ B. 

Similar to (5.45), we can write a recurrence relation which moves the black vertices by 
using the recurrence K^ 1 .K.K* = K*. Using the analogous steps as above but for the 
black vertices, we reduce the generating function X^eB ^> y K~ 1 (x, y) as a linear combination 
of K~ l (x, (0, 1/2)), K~ 1 (x, (2n, y 2 )) and additive terms which are all dependent on the white 
vertices. These additive terms are given by d s ides- We can substitute back each expression 
that we obtain for '^2 yeB 1 o v K~ 1 (x,y) into (5.57). Evaluating gives (5.12) in terms of the 
K~ l (x,y) where x is a white vertex and y is a black vertex, both on the boundary of the 
Aztec diamond. 

We have computed i^ _1 ((xi,0), (0,y 2 )) in Lemma 5.1 By symmetry, we have that 
K-\(x h 0), (0, y 2 )) = K-\(2n - x h 2n), (2n, 2n - y 2 ) 



which gives the last equation in (5.13). To compute \K 2n), (0, 2n — y 2 ))\, we can 

use the expression from \K~ 1 ((x\, 0), (0, y 2 ))\ provided we interchange between a and b. 
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The correct sign of K ((xi, 2n), (0,2ra — j/2)) is given in the third equation of (5.13) and 
is determined by computing the sign of one perfect matching of the Aztec diamond when 
removing the vertices (xi,2n) and (0,2n — 7/2)- From K {(x\, 2n), (0,2ra — 2/2)), we can find 
K~ l ({2n — xi,0), (2n, 2/2)) by symmetry. □ 



Appendix A. The Matrix N 
In this subsection, we give the matrix N. We have 



N(a, b, w, 2)1,1 = 45 5 (1 + wz) + a 2 6 3 (7 + 3w + 3z + 5to) + 2a 4 6(2 + w + 2 + to) 
N(a, 6, to, 2)1,2 = -46 7 z(l + wz) - a 2 6 5 z(5 + 5w + 52 + 5to + 2to 2 ) 

- a 4 b 3 z(2 + 6u> + 52 + 5to + 2 2 + wz 2 ) - 2a 6 bz(w + 2 + to) 
N(a, 6, io, 2)1,3 = -46 7 w(l + to) - a 2 b 5 w(5 + 5w + 5z + 5wz + 2w 2 z) 
(A.l) - a 4 6 3 u>(2 + 5w + w 2 + 6z + 5to + io 2 z) - 2a 6 bw(w + z + wz) 

N(a, 6, w, 2)1,3 = 46 9 to(1 + to) + a 2 b 7 wz(3 + 7w + 7z + 5wz + 2w 2 z + 2to 2 ) 

+ a 4 6 5 TO(3 + 9w + 2to + 9z + IOto + w 2 z + 2z 2 + to 2 + w 2 z 2 ) 

+ a 6 6 3 -ro(7 + 6w + w 2 + 62 + 7wz + w 2 z + z 2 + wz 2 ) 

+ 2a 8 bwz(2 + w + 2 + to) 



N(a, 6, if, 2)2,1 = 2a6 4 (l + w + to) + a 3 6 2 (5 + 7w + 2 + to) + 4a 5 (l + u;) 

N(a, b, w, 2)2,2 = -2ab 6 (2 + 2 + to + to 2 ) - a 3 5 4 (7 + 3w + 72 + 3to + 22 2 + 6to 2 ) 

- a 5 b 2 {4 + 2w + 7z + 7wz + 32 2 + 3to 2 ) - 4a 7 (l + w)z 

N(a, 6, w, 2)2,3 = -2a5 6 w(l + w + to) - a 3 6 4 u>(6 + 9w + w 2 + 2z + wz + w 2 z) 

- a 5 b 2 w(9 + 9w + 2to + 2 + to) - 4a 7 w(l + io) 
N(a, b, w, 2)2,4 = 2ab s w(2 + 2 + to + to 2 ) 

+ a 3 6 6 w(5 + 5u> + 62 + 3to + w 2 z + 32 2 + 6to 2 + w 2 z 2 ) 

+ a 5 b 4 w(2 + 5w + w 2 + 5z + IOto + w 2 z + 72 2 + 5to 2 + 2w 2 z 2 ) 

+ a 7 b 2 w(2w + 7z + 9wz + 2w 2 2 + 32 2 + 3to 2 ) + 4a 9 w(l + w)z 



COUPLING FUNCTIONS FOR DOMINO TILINGS OF AZTEC DIAMONDS 33 

(A.3) 

N(a, b, w, 2)3,1 = 2a6 4 (l + z + wz)+ a 3 6 2 (5 + w + 7z + wz) + 4a 5 (l + z) 
N(a, 6, w, 2)3,2 = -2ab 6 z(l + z + wz) - a 3 6 4 z(6 + 2w + 9z + wz + z 2 + wz 2 ) 

- a 5 6 2 z(9 + w + 9z + wz + 2z 2 ) - 4a 7 z(l + z) 

N(a, b, w, 2)3,3 = -2a6 6 (2 + «; + + u> 2 z) - a 3 & 4 (7 + 7u? + 2w 2 + 3z + 3wz + 6u> 2 z) 

- a 5 6 2 (4 + 7w + 3w 2 + 2z + 7wz + 3u> 2 z) - 4a 7 w(l + z) 
N(a, 6, to, 2)3,4 = 2ab s z(2 + w + wz + w 2 z) 

+ a 3 6 6 z(5 + 6w + 3w 2 + 5z + 3u>z + 6w 2 z + wz 2 + w 2 z 2 ) 

+ a 5 6 4 z(2 + 5w + 7w 2 + 5z + lOwz + 5u> 2 z + z 2 + wz 2 + 2w 2 z 2 ) 

+ a 7 b 2 z(7w + 3w 2 + 2z + 9wz + 3u> 2 z + 2wz 2 ) + 4a 9 u>z(l + z) 

(A.4) 

N(a, 6, w, 2)4,1 = a 2 6 3 (3 + io + z + wz) + 2a 4 6(5 + w + z) + 8a 6 6 _1 

N(a, 6, w, 2)4,2 = -a 2 b 5 (2 + z + wz + z 2 + wz 2 ) - a 4 6 3 (5 + w + 5z + wz + 6z 2 + 2wz 2 ) 

- 2a 6 6(2 + 5z + wz + 3z 2 ) - Sa%~ l z 

N(a, b, w, 2)4,3 = -a 2 6 5 (2 + w + to 2 + m + w 2 z) - a 4 6 3 (5 + 5w + 6u? 2 + z + u>z + 2w 2 z) 

- 2a 6 6(2 + 5w + 3w 2 + wz) - 8a 8 & -1 u; 

N(a, 6, w, 2)4,4 = a 2 6 7 (4 + 2w + 2z + 3wz + w 2 z + wz 2 + w 2 z 2 ) 

+ a 4 6 5 (7 + 7w + 2w 2 + 7z + 6wz + 3w 2 z + 2z 2 + 3wz 2 + 5w 2 z 2 ) 
+ a 6 6 3 (4 + 7w + 3w 2 + 7z + 9wz + 6w 2 z + 3z 2 + 6wz 2 + 5w 2 z 2 ) 
+ 2a 8 b(2w + 2z + 5wz + 3w 2 z + 3wz 2 ) + 8a 10 b~ 1 wz 
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